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We  demonstrate  synchronization  of  laser-induced  self-sustained  vibrations  of  radio-frequency 
micromechanical  resonators  by  applying  a  small  pilot  signal  either  as  an  inertial  drive  at  the  natural 
frequency  of  the  resonator  or  by  modulating  the  stiffness  of  the  oscillator  at  double  the  natural 
frequency.  By  sweeping  the  pilot  signal  frequency,  we  demonstrate  that  the  entrainment  zone  is 
hysteretic  and  can  be  as  wide  as  4%  of  the  natural  frequency  of  the  resonator,  400  times  the  II Q 
~  10“ 4  half-width  of  the  resonant  peak.  Possible  applications  are  discussed  based  on  the  wide  range 
of  frequency  tuning  and  the  power  gain  provided  by  the  large  amplitude  of  self-oscillations 
(controlled  by  a  small  pilot  signal).  ©  2003  American  Institute  of  Physics. 
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Frequency  entrainment  is  one  of  the  most  spectacular 
phenomena  in  nonlinear  vibrations.1,2  It  was  discovered  in 
the  17th  century  by  Huygens  who  remarked  that  two  slightly 
out  of  step  pendulum-like  clocks  became  synchronized  after 
they  were  attached  to  the  same  thin  wooden  board.  Later,  a 
similar  phenomenon  (injection  locking)  was  observed  in 
radio-frequency  (rf)  circuits3-5  and  laser  systems.6,7 

In  this  letter,  we  report  frequency  entrainment  observed 
in  a  microelectromechanical  system  (MEMS)  in  the  rf  range. 
Self-sustained  vibrations  of  mushroom-like  microfabricated 
resonators  (thin,  single-crystal,  silicon  disk  suspended  above 
the  wafer  by  a  silicon  dioxide  pillar  at  the  center)  induced  by 
a  continuous- wave  (cw)  laser  were  reported  earlier.8  Large 
amplitude  (~\/2  laser  light)  self- sustained  mechanical  oscil¬ 
lations  spontaneously  occur  at  the  fundamental  frequency  /0 
when  the  intensity  of  a  cw  laser  beam  focused  on  the  periph¬ 
ery  of  the  disk  reaches  a  critical  value  /^threshold-  The  self- 
oscillatory  behavior  is  attributed  to  the  automodulation  of  the 
laser-induced  thermal  stress  field  that  occurs  as  the  disk 
moves  through  the  interferometric  pattern  created  by  incident 
and  reflected  laser  beams.  These  spontaneous  vibrations  are 
subject  to  a  frequency  instability  A///0~  10-3  observed  on  a 
spectrum  analyzer  as  the  sporadic  motion  of  a  narrow  reso¬ 
nant  peak. 

We  can  establish  full  control  over  the  frequency  and 
phase  of  the  self-oscillations  of  the  disk  by  applying  a  small 
periodic  perturbation — a  “pilot”  signal — via  inertial  drive 
(provided  by  an  ac  voltage  applied  to  a  piezoelement)  at  a 
frequency  close  to  /0  or  by  modulating  the  effective  spring 
constant  of  the  disk  by  superimposing  a  small  rf  component 
on  the  laser  beam  intensity9  near  2  /0 .  As  the  frequency  of 
the  perturbation  is  tuned  close  to  /0  or  2  /0 ,  we  readily  ob¬ 
serve  a  distinct  point  when  the  phase  of  the  self-oscillatory 
motion  gets  locked  to  that  of  the  pilot.  On  the  polar  plot 
display  of  a  network  analyzer  (referenced  by  pilot  signal), 
the  locking  is  manifested  by  a  change  from  a  pattern  of  spi¬ 
rals  to  a  single  point.  Once  synchronized,  a  wide  range  of 
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frequency  entrainment  can  be  demonstrated  by  detuning  the 
pilot  signal  from  /0  (or  2/0)  while  preserving  the  locked 
state.  There  are  no  frequency  beats  in  the  synchronized  state: 
Only  vibrations  with  the  pilot  frequency  exist  in  the  system. 
In  our  experiment,  the  entrainment  range  A///0~4%  was 
demonstrated  for  a  MHz-range  oscillator. 

The  fabrication  process  for  rf  disklike  resonators  (disk 
diameter  30  ja m,  thickness  h  =  0.25  jam,  central  Si02  pillar 
1  jam  height,  0.6  jam  diameter)  has  been  described 
elsewhere.9  The  resulting  structure  (inset  in  Fig.  1)  was  at¬ 
tached  to  a  piezoelectric  actuator  and  placed  in  a  vacuum 
chamber  (F~10'7  Torr).  The  motion  of  the  disk  was  de¬ 
tected  interferometrically.10 

Figure  1  shows  the  amplitude  and  frequency  of  the  vi¬ 
bration  of  the  disk  as  a  function  of  the  pilot  signal  frequency. 


FIG.  1.  Frequency  and  amplitude  of  the  laser-induced  self-sustained  vibra¬ 
tions  as  a  function  of  the  pilot  signal  frequency.  Laser  power  Pcw 
=  650  julW,  perturbation  applied  as  an  inertial  drive  ( V^ezo=4  Vol). 


0003-6951 /2003/83(16)/3281/3/$20.00  3281  ©  2003  American  Institute  of  Physics 

Downloaded  31  Dec  2003  to  132.250.146.97.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://ojps.aip.org/aplo/aplcr.jsp 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

20Q2  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2003  to  00-00-2003 

4.  TITLE  AND  SUBTITLE 

Frequency  entrainment  for  micromechanical  oscillator 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory ,4555  Overlook  Avenue 

SW, Washington, DC, 20375 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

\PlSTR 

1 8 .  NUMBER  1 9a.  NAME  OF 

fTR  PACRQ  PR9Pn\TCTm  P  PPPQfYM 

/\D  L  i\/\U  1 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 

unclassified  unclassified  unclassified 

3 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


3282 


Appl.  Phys.  Lett.,  Vol.  83,  No.  16,  20  October  2003 


Zalalutdinov  et  al. 


© 

N 


2.0- 

1.5- 

1.0- 

0.5- 

0.0 


A 

A 


— A—  lockup 
A  ffree_up 
— T“  lockdown 
v  ^free  down 


v 

V 

V 


▼ 

T 
T 

! 

'  ▲  v  A 

v  *  T  £ 
v  4  J  A 
T  A 


%A 


A 

A 

A 

A 

A 

A 

A 

A 


'pilot' 


,  MHz 


1.20  1.21  1.22  1.23  1.24  1.25  1.26 


FIG.  2.  Entrainment  zone  as  a  function  of  the  pilot  signal  amplitude.  Laser 
power  Pcw=  1.7  mW,  perturbation  applied  as  an  inertial  drive.  A  disconti¬ 
nuity  in  the  width  of  the  entrainment  zone  at  VpieZo~  1-8  V  is  observed  when 
the  expanded  entrainment  zone  overlaps  the  areas  that  allow  different  types 
of  limit-cycle  oscillations. 

A  helium-neon  (HeNe)  laser  beam  at  Flaser=  650  jmW  was 
used  to  excite  the  disk  into  self-oscillation  (Fthreshold 
~  600  julW)  and  an  inertial  force  perturbation  was  applied  by 
the  piezoactuator  at  frequency  /pilot .  The  solid  line  repre¬ 
sents  data  acquired  while  linearly  increasing  /pilot .  After  es¬ 
tablishing  synchronization  at  /pilot=/lock  up  (~2kHz  below 
/o),  the  frequency  of  the  motion  follows  the  change  of  the 
pilot  frequency.  The  amplitude  of  the  mechanical  motion  in¬ 
creases  until  another  critical  point  is  reached  (/free  up,  28 
kHz  above  /0).  Synchronization  is  then  broken,  self¬ 
oscillations  at  the  natural  frequency  /0  are  restored  and  for 
/piiot>/free  up  the  pilot  signal  is  ineffective.  This  general  se¬ 
quence  is  repeated  on  a  down  sweep,  except  that  when  en¬ 
trained  at  /piiot</iock  down »  the  oscillation  amplitude  is  re¬ 
duced.  The  critical  values  of  the  pilot  frequency  /free  and  /lock 
are  different  for  up  and  down  sweeps,  demonstrating  the  hys- 
teretic  behavior  of  the  synchronization  for  our  microme¬ 
chanical  oscillator. 

The  width  of  the  entrainment  range  (/i0Ck— /free)  is  de¬ 
termined  by  the  power  of  the  pilot  signal.  Figure  2  shows  an 
entrainment  map  where  critical  points  (for  up  and  down 
sweeps)  are  plotted  as  a  function  of  the  ac  piezovoltage.  This 
plot  has  a  V-shape,  consistent  with  the  expectation  of  zero 
entrainment  range  at  zero  perturbation.  Due  to  the  hysteresis, 
two  entrainment  areas  are  present,  corresponding  to  up  and 
down  frequency  sweeps.  These  two  areas  overlap  above  the 
natural  frequency,  creating  an  “inner  V” — where  synchroni¬ 
zation  to  the  pilot  signal  is  assured.  After  establishing  phase 
lock,  it  is  preserved  for  pilot  excursions  within  the  “outer 
V” — boundaries  marked  by  /free  down  and  /free  up,  as  dem¬ 
onstrated  by  applying  frequency  and  phase  modulation  to  the 
pilot  signal. 

Subharmonic  or  superharmonic  entrainment  expected  at 
nf0  and  (1  /n)f0  (n  is  an  integer)11  could  be  more  attractive 
for  applications  because  the  pilot  signal  and  the  base  fre¬ 
quency  are  well  separated,  avoiding  parasitic  coupling.  We 
failed  to  demonstrate  entrainment  with  these  methods  using 
the  inertial  drive  as  a  pilot  signal.  However,  a  wide  range  of 
subharmonic  entrainment  was  observed  by  partial  modula¬ 
tion  of  the  laser  beam  intensit 
Downloaded  31  Dec  2003  to  13 


FIG.  3.  Entrainment  zone  as  a  function  of  the  pilot  signal  power.  Laser 
power  Pcw=  1.7  mW,  perturbation  applied  as  a  modulation  of  the  secondary 
laser  beam  (Ar+)  positioned  at  the  center  of  the  disk. 

Figure  3  shows  the  map  for  2/0  entrainment.  A  cw 
HeNe  laser,  Plaser=  1.7  mW  was  focused  on  the  periphery  of 
the  disk  to  initiate  self-oscillations.  An  Ar+  ion  laser  beam 
modulated  at  /pilot~  2  /0  illuminated  the  area  around  the  pil¬ 
lar  providing  parametric  excitation.  The  ac  component  of  the 
Ar+  laser  intensity  is  plotted  as  the  vertical  axis.  Hysteretic 
behavior  of  the  parametric  synchronization,  similar  to  that 
for  the  inertial-drive  entrainment  is  illustrated  by  the  inner  V, 
i.e.,  overlap  of  the  entrainment  regions  from  up  and  down 
pilot  frequency  sweeps.  The  outer  V,  the  range  of  frequency 
detuning  that  still  preserves  synchronization,  is  expanded  to 
70  kHz  for  50  /zW  pilot  signal  power.  Frequency  stability  of 
the  entrained  oscillations  seems  to  be  limited  by  that  of  the 
pilot  generator.  In  our  case,  A///0~10-9  was  observed. 

The  analysis  of  frequency  entrainment  in  our  mushroom¬ 
like  oscillators  can  be  obtained  from  the  differential  equa¬ 
tions  describing  self-oscillatory  motion:12 

1 

z+q(z-DT)  +  (1+CT)(z-DT+/3(z-DT)3) 

M  sin(  wpiezof ) ,  (1) 

T = A  PJ 1  +  (p  sin(  (omodt )]  -[a+y  sin2  2  n(z  ~z0)]~BT. 

(2) 

The  disklike  resonator  is  represented  in  Eq.  (1)  as  a  one 
degree  of  freedom  oscillator  with  quality  factor  Q , 
temperature-dependent  stiffness  (1  +  CT),  and  nonlinear  co¬ 
efficient  /3.  Time  is  nondimensionalized  by  the  small  ampli¬ 
tude  period  of  vibration.  T  is  the  temperature  increase  under 
the  laser  spot  and  z  is  the  deflection  normalized  by  the  laser 
wavelength  X.  The  DT  terms  represent  deflection  of  the  disk 
due  to  heating.  They  arise  from  the  slightly  bowed 
( ~  50  nm)  equilibrium  shape  of  the  disk  (a  result  of  the  re¬ 
sidual  stress  incorporated  in  the  SOI  wafer).  The  term  on  the 
right-hand  side  of  Eq.  (1)  corresponds  to  the  inertial  pilot 
signal  applied  by  the  piezoelement  at  copiezo.  Equation  (2) 
describes  the  heat  balance  for  the  mushroom.  The  incident 
laser  power  Flaser  can  be  partially  modulated  to  produce  a 
pilot  signal  at  a  frequency  wmod  (cp  =  modulation  depth) .  The 
term  in  the  second  set  of  square  brackets  describes  the  por¬ 
tion  of  the  laser  power  that  is  actually  absorbed  by  the  disk: 


y  at  /laser  =  2 /o.  A  constant  minimum  absorption  a,  and  oscillatory  compo- 
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FIG.  4.  Result  of  numerical  simulations  for  the  amplitude  of  the  laser- 
induced  self-sustained  vibrations  as  a  function  of  the  inertial  drive  pilot 
signal  frequency  (Pcw=650  julW).  The  amplitude  of  inertial  drive  M 
=  10-4,  modulation  cp  in  Eq.  (2)  has  been  taken  as  zero. 

nent  governed  by  the  product  of  y  and  sin2[27r(z— z0)]  model 
the  absorption  determined  by  the  laser  light  interference  be¬ 
tween  the  disk  and  wafer.  The  term  BT  describes  thermal 
relaxation  through  the  pillar.  Finite  element  analysis  (FEA) 
was  used  to  estimate  values  for  the  coefficients  in  Eqs.  (1) 
and  (2):  A  =  0.018  jnW~\  #  =  0.488  K-1,  C=3.53 

X  10-4  K-1,  D=  1.3X  10-5  K-1,  and  z0  =  0.06  for  a  mush¬ 
room  with  an  outer  diameter  of  40  jmm  and  a  Si02  pillar 
diameter  of  6.7  jmm.  The  absorption  properties  are  estimated 
to  be  a=. 06  and  y=  .26.  The  value  of  j8^ 0.4  is  estimated 
from  measured  frequency -amplitude  curves. 

Figure  4  illustrates  sample  results  of  numerical  integra¬ 
tion  for  the  case  of  an  inertial  pilot  signal  (to  be  compared 
with  Fig.  1):  amplitude  M=  10“ 4  corresponds  to  4  V  ac  sig¬ 
nal  applied  to  piezoactuator,  the  laser  modulation  cp  =  0.  The 
trajectory  of  the  mechanical  motion  z(t)  for  a  given  fre¬ 
quency  of  the  pilot  signal  was  calculated  at  each  step.  The 
frequency  and  amplitude  of  the  motion  was  then  determined 
from  the  calculated  z(t).  All  of  the  major  features:  Synchro¬ 
nization,  hysteretic  behavior  of  frequency  entrainment,  and 
amplitude  dependence  of  the  phase-locked  vibrations  are 
demonstrated  by  the  numerical  solutions. 

Qualitatively,  the  observed  frequency  entrainment  is 
similar  to  the  entrainment  of  a  driven  Van  der  Pol  limit  cycle 
oscillator.3,13  The  major  difference  of  the  wide-range  hyster¬ 
esis  of  the  frequency  entrainment  observed  in  our  experi¬ 
ments  can  be  explained  by  the  nonlinear  behavior  of  our  disk 
oscillator11  due  to  the  large  amplitude  motion.  Accordingly, 
the  nonlinearity  determined  shape  of  the  entrainment  zone 
can  be  tailored  by  the  cw  component  of  the  laser  power  Pcw . 
The  excess  of  Pcw  over  Threshold  enhances  the  amplitude  of 
self-oscillations  and  due  to  hardening  of  the  disk,  increases 
the  frequency  of  the  vibrations.  The  entire  entrainment  map 
(Figs.  2  and  3)  can  be  shifted  to  the  left-  and  right-hand  side 
by  adjusting  Pcw.  The  asymmetry  of  the  entrainment  zone 
(up  versus  down  frequency  detuning),  demonstrated  by  Figs. 
1  and  4,  is  defined  by  setting  Pcw  just  above  ^threshold*  where 
a  small  frequency  decrease  will  bring  the  resonator  down 
along  the  nonlinear  amplitude -frequency  curve1  toward  zero 
amplitude  and  thus  break  the  entrainment.  By  increasing 
Pcw>#threshold,  one  can  reach  a  self-limiting  regime8  when 
the  amplitude  of  self-oscillations  is  restricted  by  the  inter¬ 
ferometric  pattern  and  is  largely  decoupled  from  the  fre¬ 


quency.  This  regime  provides  a  mostly  symmetric  shape  for 
the  entrainment  zones  illustrated  by  Figs.  2  and  3. 

A  detailed  study  of  frequency  entrainment  in  MEMS  is 
in  progress  in  our  lab,  motivated  by  possible  applications  for 
high-frequency  signal  processing.  Synchronization  of  limit 
cycle  oscillators  by  parametric  excitation  (2/0  laser  modu¬ 
lation  in  our  experiment)  is  the  subject  of  recent  theoretical 
research,14,15  indicating  the  possibility  for  chaotic  behavior. 
A  theoretical  model  was  recently  proposed  for  mechanical 
neurocomputing  based  on  synchronized  MEMS  oscillators.16 
In  wireless  devices,  for  phase  or  frequency-modulated  data 
transmission,  an  rf  MEMS  oscillator,  parametrically  en¬ 
trained  by  a  pilot  signal  from  a  receiving  antenna,  would 
perform  as  a  power  amplifier,  filter,  and  frequency  down 
converter.  Synchronization  could  be  useful  for  MEMS  appli¬ 
cations  (for  example,  optical  modulators)  requiring  wide 
range  yet  extremely  precise  periodic  mechanical  motion. 

In  conclusion,  synchronization  of  a  rf  MEMS  limit  cycle 
oscillator  was  observed  by  applying  a  perturbation  signal 
either  as  a  direct  drive  or  as  a  parametric  excitation.  A  wide 
range  of  frequency  entrainment  (4%  of  the  natural  frequency 
/o)  was  observed  with  pronounced  hysteretic  behavior  attrib¬ 
uted  to  the  strong  nonlinearity  of  the  resonator.  Frequency 
stability  A///0~  10“ 9  was  demonstrated  for  entrained  oscil¬ 
lations  and  simulated  using  a  simple  model. 
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